Brazilian Journal of Probability and Statistics 

2010, Vol. 24, No. 3, 479-501 

DOI: 10.1214/09-BJPS105 

© Brazilian Statistical Association, 2010 

Reduced long-range dependence combining Poisson 
bursts with on— off sources 

David A. Rolls 

University of Melbourne 

Abstract. A workload model using the infinite source Poisson model for 
bursts is combined with the on-off model for within burst activity. Burst 
durations and on-off durations are assumed to have heavy-tailed distribu- 
tions with infinite variance and finite mean. Since the number of bursts 
is random, one can consider limiting results based on "random centering" 
of a random sum for the total workload from all sources. Convergence re- 
sults are shown to depend on the tail indices of both the on— off durations 
and the lifetimes distributions. Moreover, the results can be separated into 
cases depending on those tail indices. In one ease where all distributions 
are heavy tailed it is shown that the limiting result is Brownian motion. In 
another case, convergence to fractional Brownian motion is shown, where 
the Hurst parameter depends on the heavy-tail indices of the distribution 
of the on, off and burst durations. 

1. Introduction 

Workload models for packet traffic have been described previously [Brichet et al. 
(1996, 2000); Levy and Taqqu (2000); Mandelbrot (1969); Taqqu and Levy (1986); 
Taqqu, Willinger and Sherman (1997); Willinger et al. (1997); Kurtz (1996)]. An im- 
portant model has been the strictly alternating on-off model with heavy tailed on or 
off times [Taqqu, Willinger and Sherman (1997); Willinger et al. (1997)], in which each 
source creates work at constant rate for all time. Another important model has been 
the "M/G/oo queue" described by Cox (1984), also called the "infinite source Poisson" 
model [Mikosch et al. (2002)]. In this model, traffic arrives as independent bursts of 
heavy-tailed size or duration at the time points of a homogeneous Poisson process. There 
is no variability within each burst and bursts are considered to have the same constant 
rate. (See Mikosch et al. (2002) for a treatment of convergence results in both models.) 
The use of heavy tails is motivated by, for example, empirical evidence on the sizes of 
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WWW objects [Crovella, Taqqu and Bestavros (1998)]. For these models, the cumula- 
tive work centered about its mean is shown to be a fractional Brownian motion (IBm) in 
appropriate limiting regimes, although other limiting regimes have been studied which 
show convergence to stable Levy motion [e.g., Taqqu, Willinger and Sherman (1997); 
Mikosch et al. (2002); Kaj and Taqqu (2007)]. 

A number of more recent models combine Poisson arrivals of bursts with assumptions 
on the burst volume, burst rate, or by introducing some dynamics within the burst. In 
Maulik, Resnick and Rootzen (2002) the infinite source Poisson model is used but with 
an independent random rate given to each burst. There is no variability within bursts. In 
Maulik and Resnick (2003) bursts are given a "transmission schedule" according to an H- 
self-similar process with nondccreasing, cadlag paths (necessarily H >\), and a "volume" 
of data (work) whose distribution has infinite variance and finite mean. Under a "fast 
growth" condition they show convergence (of the finite-dimensional distributions) of the 
cumulative work, centered by its mean, to a Gaussian self-similar process that generally 
lacks stationary increments. So, their limit is generally not fractional Gaussian noise. In 
(Jaglar (2004) "flows" arrive according to a Poisson process, have infinite variance Pareto- 
distributed holding times, but packets within a flow arrive according to a compound 
Poisson process, and packet sizes have finite variance. The increment of the cumulative 
workload is shown to be fractional Gaussian noise, with Hurst parameter H depending 
on holding time. 

Cluster models have been proposed to model the packet arrival process. While not 
modeling the workload per se (absent a model for packet size), they also model variability 
within a "flow." In Hohn, Veitch and Abry (2003) and Fay et al. (2006) flows arrive 
according to a Poisson process. The number of packets in each flow and the times between 
packets arc random and may have infinite variance. Long-range dependence arises when 
the number of packets has infinite variance and finite mean. 

This paper studies the workload generated in a model where sources arrive at Poisson 
time points and have heavy-tailed "session" durations. During its session a source is 
stationary with independent on and off durations. The on times are independent and 
identically distributed (i.i.d.), as are the off times, and at least one of the two distributions 
is assumed to be heavy tailed with infinite variance but finite mean. Here "session" is 
a euphemism for a structure above that of the on-off behavior. Thus the model here 
can be considered as a hybrid of the infinite source Poisson model and the alternating 
on-off model, although the proof here owes more to the latter. The model studied here 
simplifies the one studied by Rolls (2003) without changing the main results, by assuming 
the number of sessions has achieved stationarity. 

The approach used in this paper can establish a (traditional) workload limit result 
where the variability of the workload about its mean is fractional Brownian motion un- 
der certain assumptions. Moreover, the Hurst parameter is the same as that from a 
corresponding infinite source Poisson model (i.e., from the model if there was no vari- 
ability within bursts). But the main point here is something else. The on-off behavior 
imposes additional variability on top of that from the sessions themselves. Under so-called 
"random centering" this additional variability is shown to have a corresponding Hurst 
parameter that depends on the indices of both the heavy-tailed session durations and 
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the on-off durations. The formula for the Hurst parameter is a new expression, different 
from the standard hmiting results from either the on-off or infinite source Poisson mod- 
els. This might arise if session start and end times were announced and this information 
was incorporated into the centering. Two cases arc considered, giving rise to either a 
fractional Brownian motion or a Brownian motion. It is noteworthy than under random 
centering, the limit can be Brownian motion even when both the session durations and 
the on or off durations are heavy tailed. 

The rest of this paper is organized as follows. Section 2 provides background and 
introduces necessary notation. Section 3 describes the alternating on-ofF model with ses- 
sion lifetimes. Section 4 establishes convergence results for this model. Section 5 extends 
the convergence results to weak convergence. Section 6 contains the result from several 
simulations. Section 7 provides conclusions and discusses possible future work. 



2. Background 

For the on-off processes we use the notation of Taqqu, Willinger and Sherman (1997). 
Let {W{t),t > 0} be a strictly alternating, stationary on-off process that takes the value 
during an off period and 1 during an on period. Let {Wi{t),t > 0}, i = 1,2, ... be 
mutually independent copies of {W{t),t> 0}. Define the autocovariance and mean of 
W{t) by 

r{t) = E[W^{t)] - (E[W(t)])^ and fiyv = HW(t)], 

respectively. 

The on and off times are nonnegative and independent of each other. Also, the lengths 
of the on periods are i.i.d. as are the off periods. Let fj{x),Fj{x), and fij, be the density, 
distribution function, and mean for the duration of on periods {j = I) or off periods 
(j = 2), respectively. 

The on and off times are assumed to be heavy tailed with regularly varying tails: 
Fj{x) = 1 — Fj{x) ^ x~°'^ Lj{x) as a; — > oo, 1 < aj < 2, where Lj > is a slowly varying 
function, j ~ 1,2. (In fact, one can assume only one distribution is heavy tailed. The 
results would be unchanged.) Then 

^ and a^^Y&v[W{t)~E[Wit)]]^ ^^^^ 



Finally, set amin =min{Q!i,a2} and to specify the indices, let {min, max) ~ (1,2) if 
ai < 0,2 and {min, max) = (2, 1) if ai<ot\. Under these assumptions Taqqu, Willinger 
and Sherman (1997) showed 

^ , r X ^ll\T!lUw,^y))du^TM^^t) 



482 



D. A. Rolls 



where >Cliin means convergence of the finite-dimensional distributions, Znit) is standard 
fractional Brownian motion, H = (3 — ami„)/2 and 

2 '^l^max (0 ^\ 

^lim 3 / 1 WQ _ WO V V'^'^J 

(See Taqqu, Willinger and Sherman (1997) for ai =a2-) 

Conditions so that r{t) has a tail regularly varying at infinity are given by Brichet et al. 
(2000) and Heath, Resnick and Samorodnitsky (1998). For the results here it is assumed 
that 

2 

r{u) = ^(3 - a„„„)(2 - a™„)ui-"'"-L,(u) (2.2) 

for some slowly varying function Lr{u). 

Analogous convergence results to fractional Brownian motion have been obtained for 
the infinite source Poisson model [Kurtz (1996)]. In those results it is not the number of 
sources, M, but rather the arrival rate of a homogeneous Poisson process, say A, that 
goes to infinity, followed by the time rescaling T . 



3. The model description 

In this section we define the alternating on-off model with heavy-tailed session lifetimes. 
For the on-off processes we use the same notation and assumptions as in Section 2. 
In particular, at least one of the on or off-period distributions is assumed heavy tailed 
so that 1 < amin < 2. Also, let {G{t),t > 0} be the mean-zero Gaussian process with 
autocovariance r{t) in Taqqu, Willinger and Sherman (1997) arising as 

^,^EkM(!)^ = G(t). 

For the sessions, we make assumptions on their arrivals and durations. Let {T\„,i, — oo < 
i < oo} be the arrival times of a rate A„ Poisson process on M, labeled so that 
T\^,o <0<T\„.ij where {A„, n = 1, 2, . . .} is a sequence of positive constants such that 
A„ — > oo as n — > oo. For the durations of the sessions, let {Vi} be an i.i.d. sequence with 
continuous distribution H and finite mean /iy. Here Vi will be the lifetime (i.e., holding 
time, session length) for source z. It is assumed that the distribution function H{x) of 
the lifetimes is Lipschitz continuous (e.g., satisfied if H has a bounded derivative), and 
has a regularly varying tail so we may write 

H{x) = X~°'""Lvix), asess > 0, 

where Lv{x) is slowly varying. Consideration of the integrated tail of H{x), Hi{x) = 
H{z) dz will be necessary. If Usess > 1, by Karamata's theorem [Bingham, Goldie and 
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Teugels (1987), p. 28] 



Hi{x) ^ -Lv{x) asx— >-oo. (3-1) 

^sess ^ 



In the boundary case agess — 1, with the additional assumption that Lv{t)/ 
t dt < oo then 

1 r Lv{t) 

■ at ^ oo as, X ^ oo. 



Lv{x) t 

The processes {Ta,,,;}, {Vi} and {VKi(t)} are all assumed independent. Together, the 
combination of the Poisson arrival process and the sequence of holding times defines the 
busy server process of an M/G /oo queueing system. 

A key idea in this paper is the use of random sums. That is, sums whose up- 
per limit of summation is random and obeys some convergence result of its own 
[Gnedenko and Korolev (1996)]. With random sums one must distinguish between non- 
random centering of the sum, and nonrandom centering of the summands (really a random 
centering of the sum). For the model presented here, convergence results for both kinds 
of centering will play a role and so both are discussed. 

Let 

and 

pt oo 



-oo 



where 1a{x) is 1 on the set A and otherwise. With this notation An{t) is the total 
cumulative work in [0,t) for the alternating on-ofF sources with lifetimes. Note that 
Bn{t) is the total cumulative work in [0,t) from the infinite source Poisson model whose 
"bursts" are exactly the "sessions" in A„(t) and sources arc on for the complete burst. 

The total cumulative work for the infinite source Poisson model with nonrandom cen- 
tering is easy to express 



B^{Tt)-iiv\nTt^ / 
Jo 



du. (3.2) 



One can show weak convergence of a rescaled version of this quantity to fBm with H = 
(3 — asess)/'^' under suitable assumptions. Similarly, for the alternating on-ofF model with 
lifetimes described here, with nonrandom centering, one can show weak convergence of 
a rescaled version of 

An{Tt) - fiwfJ-vKTt 
to fBm with the same Hurst parameter as for (3.2). 
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The focus of this paper is something different, namely 

A^{Tt)- ^iwBn{Tt)= / V (3.3) 

where 

X„,,(i) = -MH/]1[T,„..,T,„..+V.)W- (3-4) 

This quantity captures the variabihty from the on-off dynamics on top of the variabihty 
(at a larger time scale) from the session arrivals and departures. The goal now will be to 
establish the convergence properties of this process. 



4. Main results 

Here we present our main results, which depend on the tail indices amin and asess- The 
results can be separated into several cases by these values, as shown in Figure 1. Our first 
theorem is the limit result for case 3, while the second theorem is for case 4. In particular, 
note that although the tail indices of amin and Usess both correspond to heavy tails, the 
limiting result is Brownian motion for case 4. Cases 1 and 2 are left for future work. 



Theorem 4.1. For 1 < 



< 2 and 1 < 4 - 



<2 



C lim £ lim , ' = =aZH(t), 



(4.1) 




Figure 1. Differing cases arising from the tail indices Qmm and Osess- Case 3 leads to fBm, while 
case 4 leads to Brownian motion. The long diagonal line corresponds to 3 — Omm — Qsess = 0. 
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where F = (4 - a„dn - ctsess)/'^, crf„^ is defined by (2.1), 

(obsess ^min Obsess) (4 Q^mm ^sess) 

and {Znit)} is standard fractional Brownian motion with Hurst parameter H . 
Theorem 4.2. For 1 < Oimin < 2 and —2 < 2 — Ctmin ~ Obsess < ~lj Cksess < 2 

^X„j(u)du , 

£ lim £ lim Jo l^^=^oonA ) ^ ^,ZoAt), (4.2) 

T-i-oo n^oo \/A„T 

w/iere 

c= r{x)Hi{x)dx 
Jq 

and {Zo.5(t)} is standard Brownian motion. 



The rest of this section will prove these results in three steps. Section 4.1 presents limit 
results for the incremental process (i.e, for oo). Section 4.2 describes the asymptotic 
behavior of the variance of the incremental process, which is useful for finding the limiting 
process as r — > oo. Section 4.3 proves the main results. 

4.1. Limit results for the incremental process 

Let {Z{t)} be a Gaussian process such that for s <t: 
(i) E[Zit)] = 0, 

(ii) Var[Z(t)] = cr^/iy and (4.3) 

(iii) Cov[Z{s),Z{t)]^r{t- s) H{z)dz, 



where H{z) = 1 — H{z). It will be shown that to obtain limiting results as T — >■ oo, {Z(t)} 
takes the place of the incremental process {G{t)} in Taqqu, Willinger and Sherman 
(1997). 

Theorem 4.3. For the processes {X„,,(t)} of (3.4) and {Z{t)} of (4.3) 

C lim = Z{t). 

Proof. For any d G Z+, let (ci, . . . , Cd) be arbitrary and let < ii < ^2 < • • • < be a 
partition of the time axis. To simplify the exposition, let Ik = [tk,tk+i), k — 1,2, . . . ,d—l, 
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/o = (-00, ii), Id = [id, 00) and let 



Classify an arrival as type (i, j) , i — 0^ . . . ^d, j — i, . . . ,d, if it arrives in li and ends in 
Ij. Let Nn^ij{t) be the number of type arrivals by time t,t>0. Conditional on an 
arrival at time y, the probability an arrival is of type P(^ijj{y), is 



fo, 

H{tt+i - y), 
H{t,+i -y)- H{t, 
H{td~y), 
1, 



yih, 

y e li, i^j ^d, 
y), y^li,i<j,j¥'d, 
yell, i<j, j^d, 
yell, i=j = d. 



By Proposition 5.3 of Ross (2000, p. 273), for fixed t, the random variables Nn,i.j{t), 
i ~ 0, . . . ,d, j = i, . . . ,d are mutually independent and Poisson distributed with known 
mean (and variance) 



E[7V„,.,j(t)] = Var[7V„,,,j(i)] = A„ / P(,,,)(?/) dy = A„p(,,,)(t), 

i = Q,...,d,j = i,. 



where p(ij){t) = P(^ij){y) dy. (See RoUs (2003) for details.) For each process k 
of type it contributes to Z„ during its holding period at ti^i,...,tj. Let 

(^n,fe^(*»+i)'---'^«,fe^(*j)) be copies of (w (U+i) - ,W{tj) - mutually 

independent in i, j and k. Notice that for fixed i and j the sequence of random vectors 
is i.i.d. Then the contribution of the kth type process is 

c,+iX^f(t.+i) + ...+c,X^f (i,) 
and there are Nn^i.j(td) such processes. Thus we can write 



d-l d 
1=0 j=i+l 



j{td) I 



'An 



(Strictly speaking there is also a term from arrivals occurring exactly attk, k — 1,2,. 
Since these arrivals have probability zero that term can be safely ignored.) 
By Lemmas 8.2.2 and 8.2.3 of RoUs (2003) we have 



,d. 



k=l 



'■ — - — > z 



as n 



■ 00, 
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i = 0,...,d^l,j = i + l,...,d, where Z'*'-'^ are mutually independent random variables 
such that 



N 0,p(,-,)(id) 



2 2 



where r(u) is the autocovariance of {G{t),t > 0} and so 



d~l d d 
-P r^C- i-i D ■ 



Here {G^^'^\t),t> 0} are independent copies of {G{t),t > 0} in i and j. Since (ci, . . . ,Cd) 
and (ii, . . . jtd) arc arbitrary, the result follows from the Cramer- Wold theorem. □ 



4.2. Variance and covariance of the integrated process 

The goal of this section is to establish the asymptotic behavior of the variance V{Tt) 
Va.r[YTt] as T— >■ oo of the integrated process 



Yt= / Z{u)du. 
Jo 



(4.4) 



Since {Z{u)} is a Gaussian process with mean zero, so is {Yt}- Since {Z{u)} is stationary, 
{Yt-,t> 0} has stationary increments. 

It is helpful to relate V{t) exactly to the autocovariance r{u). This is the content of 
the following lemma. Notice in particular that since < H{s) < 1, /i is a finite constant 
for any fixed X since it docs not depend on t. 

Lemma 4.1. For the variance V(t), any constant X >0, and t>X 

V{t)=h+h + h, 

where Ii ~ 2 r{x)Hi{x) dxdy, I2 — 2cx{t — X), cx — r{x)Hi{x) dx, and I3 = 

'^IxIl^i^)Hiix)dxdy. 

Proof. Since Z(u) is a mean zero process it can be shown than 

t ry 



t nt 

Z{u) du I Z{v) dv 
Jo 



V{t) = E 

Now separate the various regions of integration. 



= 2 1 / r{x)Hi{x)dxdy. 
Jo 



□ 
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The following theorem establishes the asymptotic results for the variance V{t) as t — ?> oo 
which is needed to prove the main results. 

Lemma 4.2. Assume the autocovariance r{u) satisfies the slowly varying function con- 
dition in (2.2). For the processes {Z{t),t > 0} and {Yt,t > 0} defined in (4-3) and (4-4), 
respectively, and with 

a — 



(obsess ^min ^sess^{^ ^min ^sess^ 

Case i (l< amin < 2, -1 < 2 - amm - "sess < 1; asess > Ij- 

V{t) - <T^t^-°''-'^-°'-^"Lv{t)Lr{t) ast^oo. 

Case 4: (1< amm < 2, -2 < 2 - a,nm ~ CHgess < -1, CUsess < 2): 

lim / , , — - = lim ' 



T^ooTLv{T)Lr{T) T^oo Lv(T)Lr{T) 

V{t) ^ 2ct as i oo where c= I r{x)Hj{x) dx. 



Boundary between cases 3 and 4 < ctmin < 2, 2 — amin ~ cisess = ^1 )■ Assume c = 
r{x)H J {x) dx < oo. Then 

V{Tt) 

lim ■ — ^ — - = oo and V(t) ^ 2ct as t — > oo. 

T^ooTLv{T)Lr{T) ^' 

Proof. By Corollary 1.4.2 of the Characterization theorem [Bingham, Goldie and Teugels 
(1987), pp. 17-18] there exists X >Q such that Lv{x) and Lv{x)Lr{x) are locally 
bounded and locally integrable on [X, oo). Since only the tail of V{t) will be important, 
it is assumed t> X.ln particular, this means the integrals Ii and I2 arc nonzero. 

Case 3: By assumption 1 < agess < 2 and Ugess < 3 — amin- For /a, let Ji{y) be the 
inner integral 

ry 

Ji{v) = / r{x)Hi{x)dx 
J X 

so that I3 = 2 Ji{y) dy. By (3.1) and Karamata's theorem [Bingham, Goldic and Teugels 
(1987), p. 28], as t^c» 

T O'limi'^ ~ '^min){'2 ~ Oimin) .4-a„i„-Q 



{^sess 1)(^ ^min 0^sess^{4 ^min *^sess) 



Since 4 — amin — ctsess > 1, ^3 is regularly varying with higher order than Ii and I2 (which 
are constant and linear in t, respectively), the result follows. 
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Case 4: By Karamata's theorem [Bingham, Goldie and Teugels (1987), p. 28] 



/2 + /3 = 2 / / r{x)Hi (x) dx dy 
J X Jo 

= 2c{t-X)-2 / / r{x)Hi{x)dxdy, 



(4.5) 



IX Jy 

where c= r{x)Hj{x) dx < oo since 2 — amin ~ ctsess < — 1, as y — >■ cx3 

' rix)Hjix) dx ^ ^7"'"'^^^ -^^n^n){2 - a„.„) ^3^a,„.„-.,.„^^(^) j^^(^) 
-^y^sess -^JW ^niin ^sess) 

and so for the integral in (4.5) 

hm — I — — , . , , 

t^oo t'^-"-"''^-<^'-"Lv{t)Lr{t) 

'^i^^sess 1 )(3 ^min CKsess)(4 O^min *^sess) 

Since 4 — — asess < 1 the linear term of (4.5) is regularly varying with higher order, 
and V{t) ^ 2ct as t — > cx) while 

lim , — r = lim 



T^ocTtLv(T)Lr{T) T^oo Lv{T)Lr{T) 

The second limit requires more information about Lv{T) and Lr{T). Both positive con- 
stants and the logarithm function are examples of slowly varying functions, and they 
would give quite different limits. 

Boundary case: By assumption, 2 — amin — ctsess ~ —1 and c — r{x) x 
Hi{x) dx < CO. Since 

r{x)Hi{x) - ~rT^^ ~ - amin)x~'^ Lv{x)Lr{x), 

^ ^sess 

we know Lv{x)Lr{x) — >■ as a; — >■ oo since otherwise would make the tail too heavy for 
integrability. 

Let J2(y) = r{x)Hi{x) dx and cx = /q r{x)Hi{x) dx. Then 

V(t)=h+2 [ {c~Uy))dy, 

JX 

and so the asymptotic properties of V{t) depend, in part, on those of J2{y) as y — ^ oo. For 
J2(y), limy_j,oo J2{y) = 0, by [Bingham, Goldie and Teugels (1987), Proposition 1.5.9(b)] 
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J2{y) is slowly varying with 



,. J2{y) 

lim -— -— ■ = oo, 

■y->oo Lv(y)Lr(y) 



and by Karamata's theorem [Bingham, Goldie and Teugels (1987), p. 28] it follows that 
lim — — — - — = 1 , lim ^ = oo and 

t->oo tJ2\t) t^<x> tLv\t)Lr{t) 

hm ik^M^ . 0. 

f oo t 

The results follow immediately. □ 



4.3. Proof of main results 



Proof of Theorem 4.1. By Theorem 4.3, 

is Gaussian, mean zero and has covariance given by Gov[Z [s) , Z [t)] = r{t — s) H{z) dz, 
s <t. Now, for any t>0 and T > let Y{t) be defined by 

Yit) ^* 



^LviT)Lr{T)TH 
Then for the characteristic function 

E[exp(..r(T.))]=exp(-^^-Jg|^^ 



exp 2™- — ^ -s ] asT->cx3 

^' 2H{2H-l) ' 



smce 

T^oo LviT)Lr{T)T^" ~~ 2H{2H-1) 
by Lemma 4.2. Thus, for fixed t>0 

Y{Tt)^Yi asT^cx), 



where 



2i/(2iJ-l) ^ 
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For the covariances, by Lemma 4.2 and stationary increments 

Gov[Y{Ts),Y{Tt)] 

_ 1 iV{Ts) + V{Tt) - V{T(t - s))) 

~ 2 Lv(r)L^(r)T4-"™»-".™ 

Cr2(.s4-""""-"— + i4-a„.„-a„„ _ _ g)4-a„„.-a_. aS T CX3 

Thus, by the definition of fractional Brownian motion [Samorodnitsky and Taqqu (1994)] 
the left and right-hand sides of (4.1) have the same finite-dimensional distributions. □ 

Proof of Theorem 4.2. The proof is similar to that for Theorem 4.1 but shows, for 

^(0 = ^ 

and fixed i > 

Y{Tt) 4 y2 as T oo where Y2 - N(0, 2ct) 

and for s < t 

Cov[y(T.),y(TO] = i(nT.) + nTO-nT(t-.))) ^2^^ ^^^^^ 

so the left and right-hand sides of (4.2) have the same finite-dimensional distributions. □ 

5. Weak convergence 

In this section, weak convergence is established for both the first limit (71 — > 00) providing 
the incremental process, and the second limit (T— > 00) providing the limiting process. 

5.1. Weak convergence for the first limit 

Theorem 5.1. The convergence in Theorem 4-3 can be strengthened to weak convergence 
in the space D[0,oo) equipped with the Ji topology, and the limiting Gaussian process is 
almost surely continuous. 

Proof. Convergence of the finite-dimensional distributions was established in Theorem 
4.3 so it remains to prove tightness. For any M > 0, take < < ^2 < ^3 < M, t^—tiKl 
and define 

U^{t) = ^'-1 /y-'^' and 



(5.1) 

A = (C/„(t2) - Un{tl)f{Ur^{t3) - C/„(t2))'. 
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To prove tightness in I?[0,oo) it will be shown that there exist constants Ki > and 
K2 > (possibly depending on M) such that 

E[A] < Ki{t3 ~ h)^ , 0<<3-<i<l and (5.2) 
E[U„{t2)-Un{ti)]<K2{t3^h), 0<t3-ii<l. (5.3) 

Then the resuh wiU follow by a result due to Whitt (2002, pp. 226-227). 

First we prove the bound of (5.2). Establishing the bound of (5.3) will require little 
additional work. Let {xl^'i\ti+i), . . . , x'^'l\tj)) and Nn,ij{t) be defined as in the proof 
of Theorem 4.3 (so mutually independent in i.j, and k). For simplicity define Nn^ij = 
Nn.i,j{h)- Define 

and 

1 

Since E[Xn^u {t)] = 0, by conditioning arguments it follows that E[Ui_j^k\ = -EfTij-.^] = 0. 
In what follows below both t/ij.fe and T^^ fc will contribute to -^[A] through the nonzero 
higher moments of the summands. In addition, Ti^^k will contribute through the temporal 
correlation within each summand. 
It can be shown that 

Un{t2) ~ Unih) = 10,2,2 + To, 3, 2 + C/i,2,2 + t/l,3,2 - t/o,l,l ^ud 
Units) — Unitl) = To, 3, 3 — T'1,3,3 + [/2,3,3 ^ ^^0,2,2 ^ ^^1,2,2- 

Then the expansion of A in (5.1) has 25^ terms before simplifications. Those terms whose 
factors differ in either the first or second indices {i or j) are the product of independent 
terms. Since an individual term has expected value zero, -/^[A] gets contributions from 
only 31 nonzero terms. That is, 

i?[A] = i?[([/f,2,2 + t/f,3,2 + ro%,2 + Ul,^, + Tl2^2) 

^ (t^l,2,2 + ^0,2,2 + ^2,3,3 + ^1,3,3 + ^0,3,3)] 
+ 4i?[ — ro,2,2t^l,3,2T'i,3,3f7o,2,2 " 7o,3,2C/l,2,22o,3,3[/i,2,2 (5-4) 
~ C^l,2.2C^l,3.2ri,3,3[/i,2,2 + 7o,2,2 t^l,2,2C^,2,2 t/l,2,2 
^ To, 2, 220,3,271), 3, 3[/o, 2, 2 + '7o,3,2[/l,3,22o,3,3ri,3,3] . 

The idea now will be to bound both a"^, — r[u) and pijit^) (recall E[N„,ij] = XnPijits)) 
suitably. 
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From now K will be an arbitrary constant whose exact value will vary from one 
calculation to the next. By the assumed Lipschitz continuity of H{x) we have 



Pi,2(i3)= / H{t3-s)-H{t2~s)ds 



< / Kit3-t2)ds<K{t3-ti)^ 



Since < H{x) < 1, poAh), ^0,2(^3), Pi, 3(^3)- ^2,3(^3) are bounded above by {t^ - ti). 
Thus we have 



and for (^,J)^(1,2) 



= —E[Nn,l,2CTtv] 

<K{t^~tif 



E[Ul^^,]<Kp,,,{t,)<Kit,~t,). 



Since E[Nl,^j-N, 



n.i,j\ = ^nPi j it follows that 



1 



= ■^EiNr.^.^^EiXlS,)] + 3(7V2 1 2 - N^,i,2)aU < Kit, - t,)' 

To understand contributions from Tij-.fc notice that we can write 

2 ^u\_..rn\ Ail 



'^w - r{t) = r(0) - r{t) 



Ml +M2 



-(l-TTllW), 



(5.5) 



(5.6) 



where 7rii(t) — P{W{t) = \\W{Q) — 1), which is continuous in t since the on and off 
distributions are assumed continuous. Now by Taqqu, Willinger and Sherman (1997), 
equation (13), 



l-^ii(i) 



1 



A'l 7o 



Fi{u)du+ / Fi{t ~ u)hi2{u) du, 



where H12 is the renewal function for the interrenewal distribution Fi * F2 (i.e., H12 = 
YlkLiiEi * -^2)*^) with the density hi2{u). Using the Fundamental Theorem of Calculus 
for the first term, and Laplace transform arguments for the second term, 

^(l-^ii(t)) = — A(i)- / fi{t~u)hi2(u)du. 
dt in Jo 
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On the right side, the first term is bounded between and The second term is also 

bounded [Bhat (1972), p. 167]. Therefore, since (1 — 7rii(t)) is continuous with a bounded 
density, it is Lipschitz continuous. Using (5.6) it follows that 



\<jl,-r{u)\<Ku. 



It can also be shown that 
and by expanding the square 



-E 



- Pi.]{h)r{tk - t(fc-i)) 



<K{h-hf. 



(5.7) 



Notice that since \a^/ — ''(u)l < 2cr^- the looser bound E[Tf^ j,] < Kpi,j{t^) is also avail- 
able. 

Exploiting the idea that since {VF(t)} is a 0-1 process we know for the joint moments 
E[W''{ti)W'{tj)] = E[W{ti)W{tj)] for r, s € Z+ it can also be shown that 



■(1 - + 3M')r(t3 - ^2) + ip\,Si)cTlrr{h ~ tj) 



and 



Eirrr2 rr2 i ^1,3(^3)/, n 2\ 1, , \ , Pl^sih) 4 , 2 /, n 4 
^[^^1,3,2^^1,3,3] = ^ (1 - 2/^ )r{t3 - t2) + r Cr^ + PlAhhnj 



An 



An 



+ 2p2_3(t3)[r-(i3-i2)]^ 

Fourth moments of Ti,j^k can now be found. For example, we have that 

EiTtsJ = E[Ul,^, - 4Ul,^,Ui,3^2 + 6[/i%,3C^l,3,2 " 4[/l,3,3t/i^3,2 + Ul,^^] 

= ^^^f^ia^w - r{h - t2)) + 12pi,3(i3)a^(fT^ - r{t:, - h)) 

+ 12p?^3(fT^-r2(t3-t2)) 
2^^'^(^^)-Alt3-^2) 



< 



min„{An}' 

+ 12[pi,3(t3)fT^i^(i3 - h) +pl3{t3){2a^w)K{t3 - t2)] 
<K{t3-hf. 
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The inequalities above provide a means to bound the first 25 terms in the expan- 
sion of (5.4). For the remaining six terms we use the Cauchy-Schwarz inequahty and 
independence. For example, it can be shown that 

i?[ro,2,2C/l,3,2Ti,3,3C/o,2,2]<(i?[ro'2,2f/M,2])'/'(^^[7^'2,3f^O,2,2])'/' 

= (i?[ro%,2]i?[C/i%,2]i?[Ti%,3]i?[f/o%,2]^'/' 
2 



<K{t3~ti) 

Now every term in the expansion of (5.4) is bounded by K{t3 — ti)^ for some constant 
K > and (5.2) is shown. 

To establish (5.3) note that for — ti < 1 

EiiUnits) - U„iti)f] = E[{To,3,3 + C/i,3,3 + C/2,3,3 - Uo,l,l ~ C/o,2,l)'] 
< ^^^[(7^0,3,3 + C^l'3,3 + C^2,3,3 + ^^04,1 + '^O^,!)] 
<i^[(i3-il)' + (t3-tl)] 

<i^2(i3-ii), 

where the second line follows from independence for distinct indices and the third line 
follows from (5.5) and (5.7). Thus, weak convergence in £)[0,oo) and limiting almost 
surely continuous paths are established [Whitt (2002), pp. 126-127]. □ 



5.2. Weak convergence for the second limit 



Theorem 5.2. The convergence (as T^oo) in Theorems 4-1 o.i^'d 4-2 can be strength- 
ened to weak convergence in C[0,oo) with the uniform topology. 



Proof. Recall the integrated process 

Yt 



Z{u) du 



defined in Section 4.2 which has stationary increments. Since convergence of the finite- 
dimensional distributions is established above (Theorems 4.1 and 4.2) it remains to prove 
tightness, which we do using a moment condition [Billingsley (1968), p. 95]. 
Case 3: It will be shown that for < it < 1 and T sufficiently large 



E 



T{s+u) 



Yts? 



T^HLv{t)Lrit) 



ViTu) 



T^fiLv{t)Lr{t) 



for some constant K >Q and some constant 5 > 0. 
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Since V(t) - a'^t^" Lv{t)Lr{t) we have 



ViTu) 



T^f^Lv{t)Lr{t) 



< 



Ki{TuY" Lv{Tu)Lr{Tu) 
T^HLv{T)LriT) ■ 



Proceeding as in [Mikosch et al. (2002)] (or using the Potter bounds [Bingham, Goldie and Teugels 
(1987), p. 25]), there exists to such that for all u < 1 and Tu > to 

Lv{Tu)Lr{Tu) 
Ly{T)L^{T) 

with small 5 > Q chosen so that 2H — 25 > 1. Then 

V{Tu) .2H-S 

For < Tu < to, since V{t) ~ a'^t^^ Lv{t)Lr{t) for T sufficiently large we have 
V{Tu) ^ 2{TuY ^ 2{TuY+^ 



< 



T'^"Lv{T)Lr{T) - T^"LviT)Lr{T) " T^" LviT)LriT) ° 



t}r'u^+' 



Since 2H — 5 < 1 + (5 we have 



V{Tu) 



T^"LviT)Lr{T) 



<Ku 



for T large enough. 

Boundary case and case 4: Since {Yt,t > 0} has stationary increments it is enough to 
show for that for u > and some constant K 



E 



rpl/2 



<Ku^. 



Now, we know Yt - N{0, V(t)) and so E[Yt'^] ^ 3V^{t), while 



V{t) = 2d - 2 



Jy 



'•{x)Hi dxdy < 2ct. 



Thus for all T > and u>0 and some constant K 

4n 



E 



Ytu 

ri/2 



tightness follows, and weak convergence is established [Billingsley (1968), p. 95]. □ 
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6. Simulation results 

In this section the hmit rcsuhs obtained above are demonstrated through simulation. 
The starting point for each simulation is the "session" which arrives according to a Pois- 
son process with intensity A. In these simulations each session has a Pareto distributed 
lifetime with mean fj^sess and characteristic exponent asess, 1 < «sess < 2. To simplify 
the implementation, session lengths were truncated to integers. This step produced a list 
of start and end times for the "sessions." In a second step, these start and end times 
were converted to a timeseries {B{k), fc = 1, . . . , 25 x 10^}, the number of sessions alive at 
time k. (For these simulations, the first 50 x 10^ observations were discarded to achieve 
stationarity and the next 25 x 10^ observations were retained.) This corresponds to the 
number of busy servers in the infinite source Poisson model (i.e., the busy server process 
of an M/G/qo queue). 

On-off activity within the sessions was simulated in a third step. Starting with the 
same session start and end times, session activity within each session was simulated as 
a stationary on-off process with parameters ax, /ii, a2 and ^2- This produced a dataset 
{A{k),k= 1,...,25 X 10^} which represents the number of active (i.e., alive and on) 
sources for each k. Finally, a dataset {C(fc), fc = 1, . . . , 25 x 10^} was produced where 

C{k) = B{k) - A{k) 

and corresponds to the increments of (3.3) above. 

As shown in Table 1, several series, each corresponding to different parameter values, 
were simulated. For the datasets in series 1 through series 3, five simulations were per- 
formed in each series. For series 4, 10 simulations were performed. For all simulations, the 
parameters for the on-off activity were unchanged, with ai = a2 = 1.4 and /ii = /i2 = 100. 
The parameters for the sessions, namely asess, fJ-sess, and A are shown in Table 1. The val- 
ues of ai , a2 and asess were chosen mainly for illustrative purposes. Hurst parameters can 
be calculated for the infinite source Poisson model {Hisp = (3 — Q:sess)/2), the stationary 
on-off model [Taqqu, Willinger and Sherman (1997)] {Hon/off = (3 — min(ai, a2))/2), 
and the Hurst parameter Hi,urst/on/off described above in Theorems 4.1 or 4.2. The 
theoretical values are all shown in the table. 

Hurst parameter estimates for the simulated data were made using the "logscale dia- 
gram" [Abry and Veitch (1999)]. This technique divides the data into blocks of size 2^ 



Table 1. Parameter values for the simulations 





obsess 


f^sess 


A 


Ml. M2 




burst / on / off 


Hisp 


Hon/off 


Series 1 


1.2 


12000 


1 


100 


1.4 


0.7 


0.9 


0.8 


Series 2 


1.2 


120 


5 


100 


1.4 


0.7 


0.9 


0.8 


Series 3 


1.2 


1200 


1 


100 


1.4 


0.7 


0.9 


0.8 


Series 4 


1.8 


1200 


1 


100 


1.4 


0.5 


0.6 


0.8 
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Series 1 Series 2 



0.8 
0.75 

0.65 
0.6 



12345 12345 
Run Index Run Index 

Series 3 Series 4: = 16 to 21 

0.65 

0.6 ■ -r ■ 

0.55 ■ ■ 

T ' M 

I * 0.5 i i J- 1 1- 

i 1 0.45 ■ i 1 ■ 

0.4 ■ 
0.35 ■ 

12345 2468 10 

Run Index Run Index 

Figure 2. Estimated Hurst parameters witli approximate 95% confidence intervals for tlie sim- 
ulated data. 

and calculates the variance of the nth order Daubechies wavelet coefficients for each value 
of j. Values of j are sometimes called "octaves." Collectively, the points show how the 
data scales with increasing block length. H is estimated using weighted linear regression 
to find the slope of the graph over an interval of octaves chosen by the user. Besides 
being quick, this technique provides approximate 95% confidence intervals. 

Figure 2 shows the results from series 1-4. For each run, the central dot is the value 
of the estimated Hurst parameter, while the high and low bars show the ends of the 
approximate 95% confidence intcval. For scries 1 (top left), the estimate is pcformed over 
the interval j = 15,..., 21. Since the mean session length is 12,000, and logj 12,000 ~ 
13.6, j = 15 provides a starting octave beyond any artifacts associated with the mean 
session length. For series 2 (top right) and series 3 (bottom left), with shorter mean 
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Logscale Diagram, N=3 [ (}^j^)= (16,21), a-est = 0.0483, Q= 0.64217 ] 

18 
16 
14 
12 
10 

8 

6 

4 

2 4 6 8 10 12 14 16 18 20 
Octave j 

Figure 3. Logscale diagram for run 8 of series 4 (riglit). An ultimately horizontal line corre- 
sponds to H = 0.5, and is predicted by Theorem 4.2. 

session lengths, j = 14, ...,21 was used for the regression. In all cases the estimated 
Hurst parameter is within 0.06 of the value described in Theorem 4.1, and in all but one 
case the theoretical value is within the confidence interval. 

Figure 2 (bottom right) shows the Hurst parameter estimates for the series 4 data 
using octaves 16 to 21. The horizontal line corresponds to a value of 0.5, which is the 
Hurst parameter predicted by Theorem 4.2. The estimated values are somewhat higher 
than that predicted by theory, although the confidence intervals do cover the theoretical 
value. The logscale diagram in Figure 3 shows the difficulties in estimating the Hurst 
parameters with this data. One hopes to find an interval of linear increase on which to 
obtain the Hurst parameter estimate, but where does such an interval start? A horizontal 
line would correspond to a Hurst parameter of 0.5. A reasonable start is octave 16 for 
which the estimate would be 0.524. But, it would be equally reasonable to start at a 
larger octave, for which the estimate would be even closer to 0.5. All the data for series 
4 is similar, and so estimates using octaves 16 to 21 arc shown here. Moreover, since 
the trend in each logscale diagram is to flatness, the point estimates are likely to be 
consistently biased too large. Even so, the results from these simulations appear to agree 
with the value predicted in Theorem 4.2. In all cases the confidence interval covers the 
predicted value. 

7. Conclusions 

In this paper, the infinite source Poisson model for bursts has been combined with the on- 
oflF model for in-burst activity. Burst and on-ofF durations arc assumed heavy tailed with 
infinite variance and finite mean. By using convergence results for random centering of 
random sums, weak convergence of the workload to fractional Brownian motion is shown. 
The degree of long-range dependence is shown to depend on the tail indices of both the 
on-off durations and the lifetimes distributions. Moreover, the results can be separated 
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into cases depending on those tail indices. In one case where all distributions are heavy 
tailed it is shown that the limiting result is Brownian motion. 

The method of proof here (and in Rolls (2003)) uses iterated limits for first the 
arrival rate, and then the time scale, as in Taqqu, Willinger and Sherman (1997) and 
Brichet et al. (2000). As in those papers, one could study convergence with the limits 
taken in reverse order. This is left for future work. Moreover, the iterated limits are not 
entirely satisfactory, and a single limit as in Mikosch et al. (2002) and Gaigalas and Kaj 
(2003) would be preferable. Unfortunately, the change in details within the proof would 
make it essentially a new proof, and is left for future work. 

Random centering in traffic models appears to be a new idea. The application is not 
necessarily targeted at current TCP/IP data networks but simply any network where 
the start and end of some kind of "sessions" are signaled. It would be interesting to see 
if in a network element this can be practically used to advantage, possibly leveraging the 
"nicer" queueing behavior of Brownian motion workloads. 
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